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Abstract We present a high-angular-resolution molecular line and millimeter continuum study of the massive 
star formation site IRAS 05358+3543. Observations with the Plateau de Bure Interferometer in CO 1-0, SiO 
2-1 and H^^CO^ 1-0 reveal at least three outflows which cannot be separated in single-dish data. Observations 
at millimeter and sub-millimeter wavelengths from the IRAM 30 m telescope and the CSO provide additional 
information on the region. The most remarkable feature is a highly coUimated (coUimation factor ^ 10) and 
massive (> 10 Mq) bipolar outflow of ~ 1 pc length, which is part of a quadrupolar outflow system. The three 
observed molecular outflows forming the IRj45 05358+3543 outflow system resemble, in structure and collimation, 
those typical of low-mass star-forming regions. They might therefore, just like low-mass outflows, be explained by 
shock entrainment models of jets. We estimate a mass accretion rate of ~ 10~* Mq /yr, sufflcient to overcome the 
radiative pressure of the central object and to build up a massive star, lending further support to the hypothesis 
that massive star formation occurs similarly to low-mass star formation, only with higher accretion rates and 
energetics. In the millimeter continuum, we flnd three sources near the center of the quadrupolar outflow, each 
with a mass of 75-100 Mq. These cores are associated with a complex region of infrared reflection nebulosities 
and their embedded illuminating sources. The molecular line data show that SiO is found mostly in the outflows, 
whereas H^^CO'*' traces core-like structures, though likely with varying relative abundances. Thermal CH3OH 
comprises both features and can be disentangled into a core-tracing component at the line center, and wing 
emission following the outflows. A CO line-ratio study (using data of the J = l — 0, 2 — 1&6 — 5 transitions) 
reveals local temperature gradients. 

Key words, molecular data - stars: early type - stars: formation - protostars: individual: IRAS 05358+3543 - 
interstellar medium: jets and outflows 



1. Introduction 

It is still not known whether the physical processes leading 
to massive stars are similar to those of low-mass star for- 
mation, or whether different processes are taking place. 
Classical star formation scenarios predict m oderate ac - 
cretion! rates around lU ^ — lu " WIq yr ^ ( Shu, i!JV7 J7 
which are mcapable to overcome the radiation pressure for 



^tahler et al., 2000 ). In contrast, according to other the- 
ories following the classical accretion scenario, it is pos- 
sible to form massive stars via enhanced accretion and 



disks (lijina & Adams, 1996; 


Norberg & Maeder, 2000 


Tan & McKee, 2002 




Yorke, 2002|). For more details on 



sources more massive than approximately 10 M0 (VVolfire 



Kurtz et al. (20001 ). 



al. (2OOOD, [Richer et al. (2000D , [Churchwell (200C| ), and 



& Cassinelli, 1987). Observations over the last decade 
showed that massive star formation occurs most likely in 
a clustered mode, and theories were proposed that advo- 
cate the coalescence of protostars in dense cluster centers 



to build up the most massive sources (Bonnell et al., 1998; 
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Because molecular outflows give on large angular scales 
a wealth of information about the innermost part of 
star-forming regions, in recent years massive outflows 
have been investigated in more detail. Such outflow stud- 
ies seemed to indicate a lower degree of collimation 
for massive outflows than known for low-mass sources. 
CoUimation factors between 1 and 2 were reported as 
compared to low-mass outflows, where collimation factors 
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up to 10 are found (e.g., Richer et al. 2000 ). Based on 
those observations new outflow scenarios were suggested, 



e.g., C hurchwell (2000[) proposed that massive outflows 



maser e mission is found in the re gion with low spat ial res- 
olution (Wouterloot et al., 1988). As reported by [lofani 



et al. (1995 ) and conflrmed by ^ridharan et al. (2001 ), 
IRAS 05358+3543mm does not show any 3.6 cm emission 
down to 1 mJy. Thus, based on the total luminosity of 
the region and a regular initial mass function, it is likely 
that we are witnessing a very early stage of massive star 
formation (see ^ridharan et al. 2001 ). The spectral type 
of the most massive object of the cluster is likely between 
B2 and Bl (10 Mp < M < 13 Mp). 

' Porra s et al. (2000*), |Yao et al. (2000| ) and piang et 
|al. (2001 ) have published a series of near- infrared imag- 
ing and polarimetric studies of this region, showing that 
there are two young embedded stellar clusters, one asso- 
ciated with the 12 ^m source (age ^ 3 Myr) and one with 
IRAS 05358+3543mm (age < 2 Myr). Narrow-band ob- 
servations in the v=l-0 S(l) line of molecular hydrogen 
at 2.122 /xm delineate at least two bipolar jets emanating 
from the younger cluster, and polarimetric imaging iden- 
tifles the likely source of the outflows as being deeply em- 
bedded and undetected at ne ar-infrared wavelengths (Yac 
at al., 2000; Jiang et al., 200l|). In an accompanying paper. 



might be produced by deflection of infalling material from 
the central protostar. But the analysis of the data has 
not taken into account properly the spatial resolution of 
single-dish observations. Massive star formation sites are 
on the average far more distant (a few kpc) than well 
known low-mass sources (a few hundred pc), and thus 
their angular sizes are smaller in spite of larger linear 
sizes. Recent statistical work on massive outflows with 
11" resolution by Beuther et al. (2002b| ) shows that the 
average collimation - even with this spatial res olution - 



is high er than previously thought. Additionally, [Beuther 
et al. ( |20021^ ) flnd that low-mass correlations of outflow 



and core parameters continue up to the high-mass regime, 
suggesting that similar star formation processes are re- 
sponsible for forming stars of all masses. But their data 
also indicate that even higher angular resolution is needed 
to disentangle real source and outflow structures, and thus 
detailed studies with mm interferometers are needed. 

There have been recent interferometric studies of 
massive outflow sources , of which prominent examples 
are IRAS 201 26+4104 dCesaroni et al., 1997|, |1999D or 



G192. 16-3.82 ( |Shepherd et al., 1998^ [Shepherd fc Kurtz 



199S()7]While the outflow of IRAS 20126+4104 is rather 
collimated, G192. 16-3.82 does not show high degrees of 
collimation, and Shepherd & Kurtz (1999) propose as out- 
flow mechanism a combination of a strong, wide-angle 
wind with a weak jet. 

In this paper, we present a detailed interferometric 
and single-dish study at mm wavelengths of an extremely 
young and deeply embedded high-mass protostellar object 
IRAS 05358+3543, which is also known in the literature 
as S233IR. This s ource is part of a larger sampl e presented 



and discussed by [Sridharan et al. (20Qlj ) and [Beuther et 
al. (20l^,|,|). 



IRAS 05358+3543 has different peak positions in the 
IRAS 12 /im and 100 fim bands. While the position in 
the Point Source Catalog is based on the 12 ^m emission 
due to an infrared cluster, the 100 ^m emission peaks ap- 
proximately 40" away to the north-east and traces the 
high-mass protostellar cluster we are focusing on (Fig. 
|l|). As the mm dust continuum emission coincides with 
the 100 ^m peak, we will refer to that source as IRAS 
05358+3543mm. Based on IRAS infrared data, we es- 
timate the bolometric luminosity of the source to be 
~ 6300 L q, and the dust tempe rature to be around 47 K 
(Table | ^ridharan et al. 200lj) . The flrst CO 1-0 obser - 
vations of this region are reported by Snell et al. (1990|) , 
who detected a 15 Mq bipolar outflow and classify the 
source as a massive star formation site. The estimated 
distance of 1.8 kp c is based on obse rvations of the nearby 
Hll region S233 ( [Snell et al., 1990| ). As additional mas- 
sive star formation and outflow signposts. Class II CH3OH 
(Menten, 1991; Sridharan et al., 2001) and H2O maser 



emission (Tofani et al., 1995; Beuther et al., 2002c) is ob- 
served towards IRAS 05358+3543mm. Additionally, OH 



McCaughrean et al. (in prepj ) present new, more sensitive, 
and higher-spatial-resolution near- and mid-infrared ob- 
servations of the region, which reveal new information con- 
cerning the shocked H2 outflows and their driving sources. 

2. Observations 

We have observed /i?j4S'05358+3543mm in several tracers 
with different telescopes. Table |] summarizes the obser- 
vations described below. 

Table 1. Observation parameters of spectral line obser- 
vations; observatories are the 30 m at Pico Veleta (PV), 
the Plateau de Bure Interferometer (PdBI), the Caltech 
Submillimeter Observatory (CSO), and the Very Large 
Array (VLA); quoted system temperatures Tgyg are aver- 
age values, Av is the velocity resolution. 





freq. 


Obs 


HPBW 


Tsys 


Aw 




[GHz] 




["] 


[K] 


[km/s] 


CO 1-0 


115.27 


PV 


22 


120 


0.2 


CO 1-0 


115.27 


PdBI 


4.1 X 3.3 


300 


0.4 


SiO 2-1 


86.85 


PV 


29 


85 


3.5 


SiO 2-1 


86.85 


PdBI 


5.8 X 5.6 


150 


0.3 


hi^co+ 


86.75 


PV 


29 


85 


3.5 


hi^co+ 


86.75 


PdBI 


5.8 X 5.6 


150 


0.3 


^=^00 1-0 


110.20 


PV 


22 


120 


0.8 


CO 2-1 


230.54 


PV 


11 


250 


0.1 


CO 6-5 


691.47 


CSO 


11 


3000 


0.02 


CH3OH 


241.79 


PV 


11 


250 


0.1 



2.1. Plateau de Bure Interferometer (PdBI) 

IRAS 05358+3543mm was obse rved with the IRAM 
Plateau de Bure millimeter array ( puilloteau et al., 1992 ) 
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Figure 1. The 1.2 mm bolometer image is plotted in 
grey-scale and dotted contours (the contour levels are 
100 mJy beam~^ to 1000 mJy beam~^ in steps of 
100 mJy beam~\ in short 100(100)1000 mJy beam"^), 
and the H2 image ( [McCaughrean et al., in prep. ) overlayed 
in contours (arbitrary but linear units) . The dashed circles 
present the observed PdBI mosaic (the four dashed-dotted 
western fields were observed in only two configurations). 
The circular H2 feature at the right edge of the mosaic 
is due to the 12 fim IRAS source (black star). The white 
star shows the mm continuum peak and reference cen- 
ter of the mosaic (R.A.[J2000] 05:39:13.0 and Dec.[J2000] 
35:45:54). The offsets are given in arcsecs from the center 
of the mosaic. 



between August and October 1999 in two frequency se- 
tups]^. Five 15 m antennas equipped with dual-frequency 
receivers were used in three different configurations (4D1, 
4C2 & 4C2+N09) with baselines ranging from 24 m to 
180 m. Only four antennas were used during the August 
observations. Because of the poor summer weather condi- 
tions, the quality of the 1 mm data was not satisfactory in 
both setups, and they were only used for phase corrections. 
The 3 mm receivers were tuned to 115.27 GHz (USB) and 
86.8 GHz (LSB) to cover the CO 1-0 line, and the SiO 2-1 
and H^^CO"*" 1-0 lines simultanously. The phase noise was 
lower than 30° and atmospheric phase correction based 
on the 1.3 mm total power was applied. For continuum 
measurements, we placed two 160 MHz correlator units in 
each band. Temporal fiuctuations of amplitude and phase 
were calibrated with frequent observations of the quasars 
2145+067 and 0548+378. The amplitude scale was derived 
from measurements of MWC349 and CRL618. We esti- 
mate the final fiux density accuracy to be 15%. The 
reference center is R.A.[J2000] 05:39:13.0 and Dec.[J2000] 
35:45:54, and the wlsr is —17.6 km s~^. Ten fields (see 
Figure |^) were observed to cover the whole source, with 
the exception of the first configuration which covered only 
the 6 eastern fields. 



simultaneously, which guarantees the best alignment 
between the different maps. The on-the-fiy coverages were 
done twice with 2 sec integration time per dump and a 
4" grid (i.e., Nyquist-sampling at 241 GHz). As backends 
we used the autocorrelator and two 1 MHz filterbanks. 

Bolometer observations: 

The 1.2 mm single-dish continuum observations were con- 
ducted with the MAMBO array at the IRAM 30m tele- 
scope. F or details on the obser vations and analysis of this 
data see Beuther et al. (2002a). 



2.3. Merging the interferometric and single-dish data 

The IRAM 30 m data (CO 1-0, SiO 2-1 and Hi3C0+ 
1-0) were used to derive short-spacing information and 
thereby complement the interferometric PdBI data. The 
algorithm used to derive the visibilities corresponding to 
each pointing center of the mosaic is described by pueth 
|et al. (1996D . The single-dish and interferometer visibili- 
ties are subsequently processed together. Relative weight- 
ing has been chosen to minimize the negative side-lobes 
in the resulting dirty beam while keeping the highest an- 
gular resolution possible. Again images were produced us- 
ing natural weighting, then a CLEAN-based deconvolu- 
tion of the mosaic was performed. The final beam sizes 
are 4.7" x 3.8" (RA. 51°) for CO 1-0, and 5.9" x 5.5" 
(P.A. 65°) for H13C0+ 1-0 and SiO 2-1. The beam size 
for the 115 GHz continuum data, where only the PdBI 
data were used, is 4" x 3" (P.A. 46°). 

2.4. Caltech Submillimeter Observatory (CSO) 

We used the CSOQ 10.4 m telescope to obtain a CO 6- 
5 map of IRAS 05358+3543mm in the on-the-fiy mode 
on December, 15th 1999. At this frequency, the angular 
resolution of the CSO is 11", well suited to be compared 
with the CO 2-1 data of Pico Veleta. As backend we used 
the facility AOS. Further observation parameters are given 
in Table |l|. 



2.2. Pico Veleta (PV) 



3. Observational results 



Line observations: 

To obtain large scale information about 
IRAS 05358+3543mm, we used the IRAM 30 m 
telescope at Pico Veleta in the Sierra Nevada, Spain 
in April 1999 and mapped the whole region on-the-fiy 
in the CO 2-1, CO 1-0, ^^CO 1-0, ll^^CO+ 1-0, SiO 
2-1 and CH3OH 5o - 4o(A+) lines (near 241.79 GHz). 
Due to the continuous movement of the telescope during 
on-the-fiy observations, spurious stripes in scanning 
direction (R.A.) are found in some maps (see Fig 
Except for CO 1-0 (observed later [May 2000] to add the 
short spacings, see §^.3[) all the other lines were observed 



^ The IRAM Plateau de Bure Interferometer is supported by 
INSU/CNRS (Prance), MPG (Germany) and IGN (Spain). 



3.1. Single-dish observations 
3.1.1. Dust continuum emission 

Figure 1 presents the 1.2 mm continuum image of 
IRAS 05358+3543mm. The main peak is at the 
IRAS 100 fim peak position and has more extended emis- 
sion to the north-west and the south-west. The latter sub- 
source is directly next to the 12 ^m peak position, which is 
the center of a slightly older cluster ( Porras et al., 20do| ) 



and can be interpreted as a remnant of this older star- 
forming region. Additionally, there exists a north-western 



^ The CSO is operated by the California Institute 
of Technology under funding from the National Science 
Foundation, Grant No. AST 9980846. 
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elongation in the dust map, which will be shown t o be a s- 
sociated with molecular line emission (see section ^.1.2 ). 

Assuming that the mm continuum is mainly produced 
by optically thin dust emission w ith a grain emissivity 
index P = 2, Beuther et al. (2002a ) estimate the total gas 
mass M to ^ 600 Mq and the column density Nn^ to a 
few times 10^** cm~^ (Table ||), which shows that we are 
really deahng with a high-mass star formation site. The 
column density converts to a visual abs orption magnit ude 
Av of ~ 1000 (Av = NujOM X 10^1, iFrerking 1982| ). If 
a central star has already ignited, its free-free emission 
could be quenched by the infalling core dWalmsley, 19951 ). 



Figure 3. The ^^CO 2-1 spectrum observed with the 30 m 
at the center position of the main mm core. The system 
velocity wlsr is —17.6 km/s. 



Figure 4. The integrated CO 6-5 emission in contours 
(v=[-30,-3] km s-\ levels 90(36)270 K km s'^) is over- 
layed on the grey-scale 1.2 mm bolometer map (levels 
100(100)1000 mJy beam'^). The das hed li nes show the 
western H^'^CO"'" peak as outlined in § 3.2.3 (1 channel at 
— 13.5 km/s is chosen, levels 25(25)225 mJy beam~^). The 



overlayed raster of arrows and markers outlines the differ- 
ent outflows and sub-sources as described in Figure ||. 



Table 2. Physical parameters of the obse rved region from 
1.2 mm c ontinuum data (M an d N^^, § ^.l.l]) , IRAS {L 
and Tdust , Bridharan et al. 2001 ) an d CO 2-1 observations 
(Mout, Mout and dynamical age t, § ^.1.2| ) 



Mo 



M [Mo] 

L[L0] 
Tdust [K] 
Mout [M©] 
Mout [Me/yr] 

My] 



610 
4.0 X 10^'' 
6300 
47 
20 

6 X 10"* 
36000 



Mo 



where sizcb and sizCr are the areas of the blue and red out- 
flow lobes, respectively, mn^ the mass of the H2 molecule, 
and fmaxb and Vmax^ the maximum velocities observed in 
each line wing. A more detailed description of how the 
outflow parameters are determined is given in Beuther et 



al. (2002b). According to Cabrit &: Bertout (199C) derived 



3.1.2. Molecular emission 

Figure || presents single-dish data obtained at the 30 m 
telescope on Pico Veleta. The ^^CO 2-1 lines show 
self-absorption at the line center (Fig. and blue 
and red wing emission due to molecular outflows. Blue 
(v=[-32,-21] km s-i) and red (v=[-12,-4] km s'^) 
^^CO wing emission maps were produced by integrating 
the wing-part of the spectra. It is difficult to disentangle 
bipolar struc ture in spite of the known outflow features 
presented by Porras et al. (2000 ), but the outflow lobes 
are centered at the mm core. As outlined in ^.2.2 , the in- 
clination angle of the outflow is low, and we use the mean 
intensities of the red and blue line- wing maps in the spatial 
area known from the Plateau de Bure observations for fur- 
ther outflow parameter determination. Opacity-corrected 
H2 column densities iVb and in both outflow lobes can 
be calculated by assuming a constant ^^CO/^^CO 2 — 1 
line w i ng ratio throughou t the outflow (Cabrit & Bertout, 
1990| ). IChoi et al. (1993|) found an average ^^CO/'^CO 
2 — 1 line wing ratio around 0.1 in 7 massive star-forming 
regions, corresponding to a r(^^CO 2^1) = 0.1. We 
adopt this value for our sample as well, and we assume 
30 K as average temperature in the outflow. The outflow 
mass Mout, the dynamical timescale t, and the mass en- 
trainment rate Mout are calculated via: 

Mout = (^b X sizeb -f Nj. X sizCr) mHa 
r 

t 



masses are accurate to a factor 2 to 4, whereas the error in 
the determination of the dynamical parameters is higher, 
up to a factor 10. The derived total outflow mass Mout is 
around 20 Mq and the mass entrainment rate of molecu- 
lar gas Mout approximately 6 x 10~^ Moyr""*^ (see Table 
The latter value is derived dividing Mout by the dy- 
namical timescale t. As outlined in Beuther et al. (20d2b| ), 
such a mass entrainment rate results in accretion rate esti- 
mates Maccr around 10"'' M0yr~^. These are lower limits, 
because the outflow is not st rongly inclined to the plane 
of the sky (see section 3.2.2 ), which minimizes the wing 
emission and by that the detectable (and separable) out- 
flow emission. In spite of the low accuracy, these values 
are high and conflrm the outflow to be massive. 

All other molecular line maps presented in Figure || 
peak at positions that are clearly offset from the main mm 
core. This is in contrast with the ^^CO emission, where 
not only the wing emission is centered on the main mm 
peak, but also the integrated emission (see the ^^CO 6- 
5 map in Figure 0). All data at Pico Veleta were taken 
simultaneously, thus pointing errors cannot produce these 
offsets. The ^^CO 1-0, Hi3C0+, SiO and CH3OH images 
at the line center show a similar north-western elongation 
as the mm continuum map. Low spatial as well as spectral 
resolution of the H^^CO"^ and SiO data prevents further 
analysis of these data, but the fine wings of the CH3OH 
data show a bipolar distribution slightly shifted to the 
west with respect to the CO outflow, indicating that the 
CH3OH emission might trace a different outflow in the 
west (see §|3.2.2|). 



(Wn 



v)/2 



H. Beuther et al.: IRAS 05358+3543: Multiple outflows at the earliest stages of massive star formation 



5 



Figure 2. IRAM 30 m single-dish observations. Each image shows the 1.2 mm continuum observations in grey-scale and 
dotted contours (levels 100(100)1000 mJy beam~^) with the molecular line observations overlayed in contours (labeled 
in the top left corners). We show the integrated intensity of ^^CO 1-0 (v=[-22, -9] km s"\ levels 32(4)52 K km s~^), 
Hi3co+ 1-0 (v=[-26,-10] km s'^, levels 0.4(0.6)5.2 K km s'^) and SiO 2-1 (v=[-26,-10] km s'^ levels 
0.6(0.6)5.6 K km s~^). The ^^CO 2-1 map presents blue (solid contours, v=[-32, -21] km s'^, levels 13(9)40 K km s'^) 
and red (dashed contours, v=[— 12, —4] km s^^, levels 13(9)40 K km s~^) outflow emission. For thermal CH3OH 5o — 
4o(A+) emission at 241.79 GHz, we present one map for velocities near the Hue center (v=[— 17, —16] km s~^, levels 
1.1(0.5)4.6 K km s~^) and one map for the Hue wings (blue: solid lines v=[— 20, —23] km s"^, levels 2.1(0.5)4.6 K km s~^; 
red: dashed Hues v=[— 14, — 12] km s~^, levels 1.1(0.5)1.6 K km s~^). The stretches in east-west direction, visible in 
the ^^CO and ^^CO maps, are artefacts of the on-the-fly observing mode. The beam size for each transition is shown 
in the bottom-right corners. 



Figure 5. Presented are in grey-scale the ratio maps of 
i^CO 6-5/2-1 (left) and ^^CO 1-0/2-1 (right) 
separated in b lue (l ower) and red (upper) wing emission 
(as defined in § 3.1.2 ). In the blank white regions the signal 
to noise ratio is not sufficient to derive proper ratios. The 
contours show the 1.2 mm dust continuum map (levels 
150(150)1050 mJy beam~^). The asterisk marks the hot- 
blue position (PI), the triangle the hot-red position (P2), 
and the circles the coldest position (P3). 



3.1.3. Temperature gradients 

Multi-line studies of rotationally excited CO provide a 
good tool to study temperature variations in molecu- 
lar clouds. In thermaHzed gas, where the local thermo- 
dynamical equihbrium (LTE) approximation appHes, the 
low-J CO 1-0 and 2-1 lines trace cooler material while 
the mid-J C O 6-5 line is sensitive to warmer molecu- 
lar gas (e.g., Beuther et al. 2000 , Hirano fc Taniguchi 
2001|)7|To get an idea of the temperature distribution in 
IRAS 05358+3543mm, we derived ^^CO 6 - 5/2 - I and 
^^CO 1 — 0/2 — 1 ra tio ma ps for the blue and red wing emis- 



sion as defined in i3.1.2| (Figure |). All data are smoothed 
to the resolution of the CO 1-0 transition (22"), even the 
6-5/2-1 ratio map to increase its signal-to-noise ratio. 

In most regions rather uniform ratios around 0.5 are 
observed, but each of the ratio maps has a prominent re- 
gion with line ratios larger than 1. The blue 6-5/2-1 ra- 
tios rise slightly east of the /Kj4 5 Position (named PI and 
marked by an asterisk in Figure ||), and the red 6-5/2-1 
ratios rise south-east of the main mm peak (P2, triangle 
in Figure P), which is not prominent in the blue wings. 
Contrasting to the increases in the 6-5/2-1 ratios around 
PI and P2, the 1-0/2-1 ratios are rather uniform around 
0.5 throughout those regions, but they rise - independent 
of the line wings - around 30" west of the mm peak, where 
the mm emission drops significantly (P3, circle in Figure 
I). The blue 6-5/2-1 ratios are lowest there. 

As the average ^^CO/^^CO 2 - 1 0.1 line wing ratio 
(^ 3.1.2 ) corresponds to a ^^CO 2-1 opacity of « 0.1, the 
average ^^CO 2 — 1 opacity in the wings is about 6 (Langer 



& Pen: das, 1990). To quantify the temperatures in those 
regions, we ran several Large Velocity Gradient models 
(LVG) for optically thick CO emission {t{^^CO 2-1) « 6). 



Figure 6. CO ratios versus kinetic temperatures calcu- 
lated via an LVG model for optical thick CO emission 
(t(12C0 2 - 1) « 6). The left panel shows the ^^CO 6 - 
5/2 — 1 ratio at densities of 10^ cm~'^, and the right panel 
presents the ^^CO 1—0/2—1 ratio at densities of 10^ cm~^. 



The regions around PI and P2 have average d ensities of 



a few times 10° cm ^ ([Beuther et al., 2002a| ), which is 



sufficient to thermalize even the J = 6 — 5 transition. To 
thermalize the CO 2-1 transition, densities of 10'* cm"'^ 
are sufficient, which is a reasonable assumption at posi- 
tion P3. Fig. ^ presents the resulting fine ratios versus 
the kinetic temperatures for the above estimated densi- 
ties, as calculated via the LVG approximation. Obviously, 
increasing ^^CO 1 — 0/2 — 1 ratios do indicate decreasing 
temperatures, whereas increasing *^C0 6 — 5/2 — 1 ratios 
are signposting enhanced temperatures in such regions. 
The LVG calculations result in temperatures > 80 K at 
PI and P2, and temperatures below 20 K around P3. 



3.2. Interferometric observations 

3.2.1. The main core at high resolution 

Zooming in the main single-dish mm core with the PdBI at 
2.6 mm at a spatial resolution of 4" x 3" (7000x5000 AU), 
it resolves into three massive sub-sources with separa- 
tions between 4" and 6". We label those three sources 
mml, mm2 and mm3 (see Figure The major and mi- 
nor core sizes were determined by fitting two-dimensional 
Gaussian to the data, the integrated emission is derived 
within approximately these areas. Assuming optically thin 
dust emission, masses and column densities presented in 
Table || were cal culated at 30 K (for m ore details on the 
calculations see Beuther et al. 2002a ). Source mml co- 
incides within less than 1" with the main mid-infrared 



sourc e (at 11.7 /im) presented by [McCaughrean et al. (in 



prep 



as well as w ith t he deeply embedd ed source found 



by |Yao et al. (2000|) and [Jiang et al. (200l|) in their K-band 
polarimetric imaging study. 



Three H2O maser features found by Tofani et al. (1995) 
are associated with mml, while one is associated with 



mm2. In 1999, Beuther et al. (2002c) detected only one 
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Table 3. Parameters derived for the PdBI mm cores: po- 
sitions, major and minor core sizes, peak and integrated 
intensities, and masses and column densities. 



mml 



mm 2 



mmS 



R.A. [J2000.0] 
Dec. [J2000.0] 



maj. Xj mm 
peak 



5:39:13.08 
35:45:50.5 
6.6 X 4.5 



5:39:12.78 
35:45:50.6 
6.6 X 4.1 



5:39:12.49 
35:45:55.2 
6.1 X 4.1 



a collimation factor of 10. A number of emission peaks 
along the outflow show red and blue emission in CO and 
SiO (Fig. ||a,b), which is a typical feature for expanding 
bow shocks near the plane of the sky. Thus, we conclude 
that the outflow is not stro ngly inclined to the plane of 
the sky. Porras et al. (2000|) claim th at IR93 is powering 



the outflow, but Yao et al. (2000|) and Jiang et al 



(200l[) show that IR93 is highly polarized and thus not a 



int. [mjy] 
M [Mo] 



30 
100 
5.1 



22 
73 
3.5 



23 
77 
3.5 



separate source. Our mm data, the mid-infrared images 
by McCaughrean et al. (in prep.), and the polarimetric 
images by [Yao et al. (200C| ) and [Jiang et al. (2001D 



Figure 7. The grey-scale shows the 3 mm 
continuum emission (mml, mm2, mm3; lev- 
els 75(25)200 mJy beam"!) with the CO 1- 
high-velocity outflow (B) overlayed in con- 
tours (blue: solid lines v=[-43,-29] km s~^, levels 
0.5(1)3.5 Jy beam~^ km s~^; red: dashed lines v=[- 
9,0] km s-\ levels 0.8(0.8)4 Jy beam^^ km s^^). 



strongly suggest that a protostar within mml 
powering source of outflow {A) . 



the 



of the previously known H2O maser feature south of mml 
(Fig. 

The CH3OH maser feature is located near the mm 
cores as we ll, but the pos ition is based on single-dish ob- 
servations ( Menten, 1991 ) and not accurate enough for a 
closer interpretation. 

Poi rel^T^aT (2000 j tound 3 mtrared sources (1K56, 



A high-velocity CO outflow (B) 

A second bipolar outflow is seen in CO at high veloc- 
ities relative to the system velocity (Fig. ^ blue lobe 
v=[-43,-29] km s'^ , red lobe v=[-9,0] km s'^). As Fig. 
shows, this outflow is driven most likely by a source 
within the same protostellar condensation as outflow {A), 
and both outflows together form a quadr upolar system, 
inclined by a position angle PA of ~ 40°. McCaughrean 



et al. (in prep. ) differently speculate that the second 
mid-infrared source 2" east might be the powering source. 
Blue and red lobes are clearly resolved by the PdBI with 
the red one to the south-east and the blue one to the 
north-west. The high- velocity outflow corresponds to 



the second outflow detected in H2 emission by Porras 
et al. (2000[) and further discussed by McCaughrean 



and min2. [Yao et al. (200(j| ) and piang et al. (200l[ ) then 
showed that the emission from IR58 and IR93 is highly 
polarized, and thus both are not independent sources but 
reflection nebulae powered by mml. 

3.2.2. Three (at least) bipolar outflows 

The most striking features of IRAS 05358+3543mm are 
the three outflows observed with the Plateau de Bure 
Interferometer in CO and SiO (Fig. ||). 

A large scale highly coUimated CO outflow (A) 

Figure ^ shows a highly collimated large scale (« 1 pc) 
molecular outflow observed in CO in the east of the 
region. We present the PdBI data without merging the 
single-dish data in Figure because the spatial flltering 
properties of the interferometer help to isolate the flow 
fro m the ambient gas (the combined data are presented 
in § 3.2.4 ). To highlight the differences between outflow 
{A) and the high- velocity outflow (B), which is presented 
in the next paragraph, we used slightly difi^erent velocity 
windows compared to the single-dish image (Fig. ^. 
Outflow (A) is slightly bent and oriented from north to 
south terminating in H2 bow shocks. This corresponds 
to one of the two outflows discussed by Porras et al. 



IR58 and 1K93 m f igure [/phat are all ottset troni mml gt al. (in prepj). It has a collimation factor around 6 



(200c ) and by McCaughrean et al. (in prep,) based on 



H2 data. To quantify the degree of collimation, we divide 
the length of the outflow by its width, which results in 



(estimated from combined CO and H2 data). At the 
tip of the south-eastern lobe, shocked SiO and H2 
emission is observed (Figure ||a,b). This is also the 
position of the warm region P2 (see section 3.1.3[ ) , which 
suggests that the temperature increase may be caused by 
shock interaction of the outflow with the ambient medium. 

A third outflow mainly observed in SiO (C) 

Further to the west, a third outflow is observed mainly 
in SiO 2-1 on large scales in north-south direction (Fig. 
^d) . Red and blue lobes are resolved with three symmetric 
bullet-like features in each lobe. But as the region is very 
complicated, we present two alternate scenarios for the 
observed outflow features in this region. 

(1) Outflow (Cbont): The flrst hypothesis of outflow 
structure is outlined by the two western arrows in Fig. |[ 
This way, the outflow shows the same bending as the large 
scale CO outflow (A) in the east, and it follows to a large 
degree the mm emission as outlined in Fig. ^. CO emission 
is found there with the same bending structure (Figure 
^). The collimation factor of this outflow is around 3, and 
the powering source is possibly the western H^'^C0+[3] 
core ( |3.2.3[ ). 

(2) Outflow (Ccavity): The other possible interpre- 
tation of the observation is sketched by the two ellipses 
presented in Fig. ||. This way, the northern SiO emission 
is not tracing the main part of the outflow, but rather the 
western cavity wall. The northern elongated CO feature 
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(Figure |^), which is not accounted for in the previous 
scenario, could then represent the bow shock at the tip of 
the cavity. A problem with this interpretation is that the 
SiO emission in the north is mainly red-shifted with regard 
to the system velocity, whereas the northern CO emission 
is on the blue side of the spectrum with a difference of 
approximately 10 km s~^. A possible solution of this dis- 
crepancy is that the outflow is in the plane of the sky at 
the front side of the dense gas traced by the 1.2 mm dust 
continuum emission. Thus, the SiO emission is produced 
at the backside of the cavity, where the outflow interacts 
with the dense gas, and the material moves away from the 
observer. The blue CO feature then could be produced in 
the final expanding bow shock, when the outfiow inter- 
acts with less dense gas at the tip of the fiow, which is in 
front of the outfiow and thus pushed towards the observer, 
producing blue-shifted emission. In this interpretation, not 
Hi3CO+[3] b ut rat her the B.^^CO+ core [2] is the center of 
the outfiow (§ 3.2.3| ), and the bending structure, especially 
of the northern lobe, is less prominent. The collimation 
factor of this structure is approximately 3 as well. 

The southern part of the outfiow seems to be less 
complicated and can be interpreted in both scenarios. 
At the edge of the more evolved IRAS 12 fim cluster, 
we fi nd the hot blue region PI described in section 
3.1.3 . The temperature increase of the outfiowing gas 
in this region could be caused by UV heating of the 
adjacent more evolved cluster. We point out that the 
details of the morphological interpretation of out- 
fiow (C) have no bearing on the overall energetics of the 
whole outfiow system, which is dominated by outfiow (A). 



Additional features 

Are there even more outfiows? We tentatively identify one 
more outfiow oriented in east-west direction. Molecular 
line emission of CO and SiO shows extensions in that di- 
rection, in the interferometric data as well as with the 
single-dish observations (dotted east-west line in Figures 
^, H, and Furthermore, at the western end of these 
molecular extensions, H2 emission is detected. One of the 
mm sources mm2 or mm3 might be the center of this ten- 
tative fiow. Finally, we note that even more outfiows are 
indicated in the H2 data of McCaughrean et al. (in prep.). 



3.2.3. emission 

The PdBI data of H^^CO+ 1-0 show three peaks in east- 
west orientation, and more extended emission to the north 
and south. One Hi3C0+ peak (Hi3C0+[l]) is associated 
with the three mm continuum peaks (mml to mm3). The 
other H^^CO"*" peaks [2] and [3] are both close to the 
center of outfiow (C), but we cannot determine with cer- 
t ainty, which of the two is the powering source of this fiow 
(§ 3.2.2 ). Assuming local therm odynamic equihbrium , an 
HCO+ abundance of 1 x 10^^ (van Dishoeck, 1993 ) and 
a C to ^^C ratio of 67 ( [Langer fc Penzias, 199C ), we can 
estimate the approximate masses of the H^'^CO"'" clumps, 
which are Hsted in Table |[ The temperatures decre ase to 
the west according to the line ratio results found in §3.1.3. 



Table 5. Masses of the H^'^CO"'" clumps at the given tem 
peratures. The temperatures decrease to the west ( ^.1.3 ) 



Outflow parameters 

We calculated the outfiow parameters from the merged 
CO Plateau de Bure and 30 m data with the same 
assumptions and within the same velocity range as 
outHned in §3.1.2. The derived values (Table ||) agree 



H"CO^ 



[1] [2] [3] 



reasonably well with the single-dish results (Table g), 
and make us confident that the orders of magnitude 
are correct. For this calculations we cannot disentangle 
properly the contributions of the outfiows (A) and {B), 
thus the derived values include both outfiows. But from 
the spatial extent of both outfiows, it becomes clear that 
most of the outfiowing mass and of the mass entrainment 
rate is due to the large scale fiow (.4). The value of 
outfiow (C) is calculated for the morphologically bent 
interpretation (Cbont)- 



Table 4. Outfiow parameters from the merged PdBI and 
Pico Veleta observations. 



Mout 

[Mb] 



[IQ-^Meyr 



Mout 
In 



t 



[yr] 



{A)+{I3) 

(Cbont) 

All emission 



9.6 
4.4 
16.9 



2.6 37000 
1.4 31000 
4.6 37000 



mass [Mq] 
T [K] 



19 
30 



6 

20 



15 



The three ll^^CO+ clumps are aligned in east- west di- 
rection with similar fiuxes in each clump. Morphologically, 
this is different from the mm continuum emission which 
decreases towards the west (see Fig. 0). Therefore, the 
H^^CO+ peaks correspond only weakly to column density 
and mass concentrations traced by the mm continuum. 
As we do not have other high-resolution molecular line 
observations, it is difficult to discriminate between abun- 
dance and density variations, and it is probable that the 
H^^CO+ peaks [2] and [3] are caused by an interplay of 
different processes. The mass of H^'^C0'^[1] is about an or- 
der of magnitude below the value we derive from the mm 
continuum (Tables | &|). Thus, in this region H"CO+ is 
likely depleted. Remarkably, we detect no mm source near 
H^^C0"''[2] and [3] and therefore at the center of outfiow 
(C). The 2.6 mm PdBI data show at the Hi3CO+[3] po- 
sition just a peak at the 2a level, and our Scr mm sensi- 
tivity corresponds to a mass limit of ~ 50 Mq (at 15 K). 
Even by a factor 5 enhanced H^^CO^ abundances result 
in lower mass estimates based on the H^'^CO"'' data, and 
the clump stays undetected in the mm continuum with 
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Figure 8. Presented are the PdBI observations as contour overlays on the grey-scale H2 data (McCaughrean et 



al., in prep. ), (a) Outflow (A): CO 1-0, red wing emission (v=[-14,-7] km s ^, levels 5(10)95% from the peak 



intensity, PdB data only), and blue wing emission (v=[-30,-21] km s~^, levels 10(10)90% from the peak intensity, 
PdB data only), the box marks the region shown in Figure |^. (b) Outflow (C): SiO 2-1, red wing emission (v = 
[—15.5,-8.5] km s~^, levels 0.15(0.15)1.15 Jy beam~^ km s^^), and blue wing emission (v=[-28.5,-20.5] km s~^, 
levels 0.15(0.15)1.75 Jy beam^^ km s^^; PdBI data with short spacings from PV). (c) R^^CO-^ 1-0, the integrated 
intensities are shown in green contours (v=[-20.5,-10.5] km s~^, levels 0.18(0.18)1.08 Jy beam~^ km s~^, PdBI data 
with short spacings from PV). The numbers in brackets label the three H^'^CO'*" sources, the beams are shown each 
time at the bottom-right. In all images, the arrows outline the three main outflows, (-4): long north-south arrows in the 
east, (B): short arrows in the east, and (Cbont): long arrows in the west. The ellipses outline the second interpretation 
of the western outflow (Ccavity), and the dotted east-west line shows the speculative forth outflow. The three squares 
represent the mm sources, the diamonds locate the H^'^CO"'" peaks [2] and [3], and the triangle marks the 12 fj,m 
IRAS position. 



our sensitivity limit. This is also the region P3, where the 
temperature decrea ses sig niflcantly according to the CO 
1 — 0/2 — 1 ratios (§ 3.1.3 ), and where only weak emission 
is observed in the integrated CO 6 — 5 map (Fig. U). 



3.2.4. Channel maps of the CO data 

In Figure ^ we show for clarity the CO 1-0 data from the 
interferometer only, since the large-scale flow is more pro- 
nounced than in the merged data. Figure § now presents 
a channel map of the Plateau de Bure data merged with 
short spacings obtained at Pico Veleta. The eastern large 
scale outflow is still visible (e.g., channels at —25.5, —23.5, 
—7.5 and —5.5 [km s~^]), and new features appear by 
adding the short spacings due to bright extended emis- 
sion at velocities close to the system velocity. 

The ring-like structures seen in a number of channels 
are no artefacts caused by the merging with the interfero- 
metric data, because the ring is also seen in the ^^CO 6-5 
single-dish data (Fig. ^). At —21.5 km the ring-like 
structure seems to surround the position of the H^^C0"*"[3] 
peak, which suggests that self absorption could cause the 
ring-like structure. But as the center of this structure 
moves spatially between the different channels, the ring 
might be produced by different overlapping outflows. The 
east- west structure at —23.5 km s~^ has the same orienta- 
tion as the putative fourth flow, which supports our idea 
of another outflow. 



3.3. Interpretation of single-dish maps with PdBI data 

Som e of t he puzzHng single-dish features outHned in sec- 
tion 



3.1.2 can be explained with the high-resolution data. 



Being near the plane of the sky, outflow (A) is difficult 
to resolve from the ambient gas in the single-dish data, 
and the outflows are far better visible with the Plateau de 
Bure Interferometer alone. Its higher spatial resolution is 
capable of resolving different flows spatially, whereas the 
interferometric feature of flltering out large scale uniform 
emission makes outflow (A) a very prominent structure 
(Fig.|a). 



Interestingly, the offsets between the molecular line 
peaks of Hi3C0+, SiO, and CH3OH, and the 1.2 mm con- 
tinuum peak are based on different physical processes in 
spite of their similar spatial distribution in the single-dish 
observation (Fig. |). In the case of Hi3C0+, the PdBI 
data resolve the single-dish map into three sub-sources 
with an extension along the north western 1.2 mm con- 
tinuum ridge, which makes H^'^CO"'" being at least partly 
a tracer of core-str uctur e, though likely with varying rela- 
tive abundances (§ 3.2.S ). In contrast, SiO is only found in 
the outflows, but because IRAS 05358+3543mm exhibits 
at least three outflows, the single-dish map looks very sim- 
ilar to the H^^CO+ image. 

The situation is different with regard to the thermal 
CH3OH emission at 241 GHz. While the line center emis- 
sion corresponds to the dust emission, the wing emission 
seems to be associated with outflow (C). From these data 
it appears that CH3OH is a tracer of outflows as well as 
of core emission. 



4. Discussion 

4.1. The 1 pc, highly collimated outflow (A) 

The eastern large scale outflow {A) in 
IRAS 05358+3543mm is the flrst example of a highly 
collimated, jet-like, bipolar and massive outflow with an 
extension of « 1 pc. The collimation factor 10 is as high 
as the highest found in low-mass flows. High degrees of 
collimation are difficult to explain by stellar winds or 
deflection of infaUing matter. Outflow models involving 
highly coUimated jets e ntraining the surroun di ng material 
are m uch more Hkely ( Cabrit et al., 1997 ). Churchwel] 
(20001) argues that the predictions of shock-entrainment 



models for the amount of material that can be entrained 
are far too low to account for several tens of solar 
masses found in massive outflows, as it is the case for 
IRAS 05358+3543mm (Table |). But he also stresses the 
caveats in these estimates, because they are highly sensi- 
ble to different assumptions, especially the entrainment 
efficiency in the dense interstellar medium as found in 
massive star formation sites. It is well possible that the 
entrainment efficiency rises significantly in the very dense 
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Figure 9. Channel map of the region in CO 1-0, Plateau de Bure data merged with short spacings from Pico Veleta. 
The grey-scale levels are between and 10 Jy beam~^ km s~^, and the contours range between 0.1 and 5.1 (by 1) 
Jy beam"^ km s~^. The bottom-right corner shows the centroid velocity of each channel. The overlayed arrows and 



markers in two channels are the same as in Fig. 3.2.2 



interstellar medium, and higher efficiencies could easily 
account for the observed high outflow masses found in 
massive star-forming regions. Therefore, outflow (A) links 
massive outflow phenomena to scenarios already known 
from low-mass star formation. It strongly suggests that 
entrainment by bow shocks propagating in a collimated 
jet also has an important role in the formation of massive 
outflows. 



4.2. The high-velocity outflow (B) and the quadrupolar 
structure 

Emanating at a PA of ~ 40°, most likely from the same 
mm core mml (Fig. the less extended high- velocity 
outflow (B) is clearly distinct from the large scale out- 
flow (A). Different quadrupolar outflow mechanisms are 
discussed in the literature (for a brief compilation and 
discussion see 



Gueth et al. 2001 ), and the most appealing 
explanation in the case of IRAS 05358+3543mm seems 
to be that the two outflows are produced independently 
by adjacent protostars inside the same mm condensation 
mml. Being parts of one and the same outflow is unlikely, 
mainly because the large scale outflow does not show any 
elongation at a PA of 40°, which would be expected if 
the quadrupolar structure was due to deflected or pre- 
cessing material. Assuming half of the HPBW (~ 1.75" « 
3100 AU) at a distance of 1.8 kpc as the the maximal sep- 
aration of the two powering sources, this is well within the 
range of typical binary separation. 



4.3. The SiO outflow (C) 

In the bent interpretation, the SiO outflow (Cbcnt) closely 
follows the north-south elongation of the dust continuum 
emission. It is slightly less collimated than the large CO 
outflow (A), and blue and red lobes are separated in SiO 
as well as in CO. The northern and southern lobes show 
three peaks each that are symmetric with respect to the 
assumed origin, suggesting three different ejection events. 
In contrast to this, the cavity interpretation of outflow 
(Ccavity) does uot imply that the outflow has to follow the 
core structure, but rather that SiO is only observed in 
that part of the flow where the densities are high enough 
to excite SiO sufficiently. As already mentioned, outflow 
(Ccavity) then is likely to be located at the front side of the 
core outlined in Fig. |l|. The bullet-like structures and their 
interpretation as three ejection events are independent of 
the exact scenarios. 

It is a rather intriguing fact that SiO, which forms 
in shocks from dust disruption (Schilke et al., 1997), is 



Q). Outflow (C) follows to a large degree the north-south 
dust fllament, while the eastern large scale CO outflow 
leaves the dense core very soon. Thus, throughout most 
parts of outflow {A), column densities are low and the 
gas densities there are not high enough to excite the SiO 
sufficiently. Another reason may be that that C-shocks 
(shock ve locities Vs < 50 k m s~^) are needed to pro- 



duce SiO (ISchilke et al., 19971) . In the less massive outflow 
(C) C-shocks are likely to be more common than in the 
more massive outflow (A), whe re higher jet velocitie s are 



expected (> 500 km s , e.g., Eisloffel et al. 200C) and 



only found in regions with strong dust emission (Figure 



shocks are more likely of J-type {vg > 50 km s~-^), disso- 
ciating the molecular material. These different scenarios 
have to be checked for a larger sample. 

4.4. The large scale bending 

Assuming the bent interpretation of outflow (Cbcnt), both 
large scale outflows and the mm dust core show the same 
bent morphology. For the cavity interpretation of outflow 
(Ccavity), Still outflow (A), the southcm lobe of outflow 
(Ccavity), and the mm dust emission exhibit the same bent 
structure. Possible bending mechanisms for protostellar 
jets are discussed, e.g., in Fendt & Zinnecker (1998). 

Internal bending scenarios on small scales seem to be 
ruled out in the case of /i?j45 05358+3543mm, because it 
is not likely that in both large scale outflows exactly the 
same small scale physics takes place (e.g., acceleration of 
the jet source by a binary component, or precession of an 
accretion disk in a binary system), and additionally the 
dust core follows a similar morphology. 

Considering the whole source structure, external ef- 
fects seem to be more plausible to explain the structure 
of this massive star-forming site. On larger scales (4 to 
10 pc) other Hii regions are found: in the north-east S233 
(at a distance of ~ 4 pc), in the west S235 (at a distance 
of « 10 pc), in the north-west S231 (distance « 4 pc) and 
a bit further away S232 (distance ~ 15 pc, for large scale 
images see, e.g., Porras et al. 200(]| ). As a whole, the region 
contains different massive star formation sites at the edge 
of IRAS 05358+3543mm, and the most likely explanation 
of the overall bending of this youngest site is due to en- 
ergy input (UV radiation as well as stellar winds) from the 
more evolved star forming regions in the vicinity. Thus, we 
are possibly observing an example of sequential star for- 
mation and interactions of regions of different ages. 

4.5. Outflows from young stellar objects of all masses 

A key question triggering this observational study was 
whether massive outflows are generated by different phys- 
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ical mechanism than their low-mass counterparts. We 
therefore compare IRAS 05358+3543mm with other out- 
flow sources of different core masses. 

A low-mass object: a prototype of a low-mass outflow is 
HH211 emanating from a dust condensation of ~ 0.2 M© 



(Gueth & Guilloteau, 1999). It includes a highly colli- 



mated, high-velocity molecular jet surrounded by a less 
collimated cavity-like outflow at lower velocity. A shock- 
entrainment model nicely explains the outflow structure. 



An intermediate-mass object: Gueth et al. (2001) re 



three massive mm cores (between 100 Mq and 75 Mq) at 
the center of the quadrupolar outflow. 

The data suggest that the physical processes associ- 
ated with this massive outflow are similar to those driving 
their low-mass counterparts. It is Hkely that many massive 
star formation sites can be shown by interferometric high- 
resolution observations to be composed of basic features 
known from low-mass star formation. The accretion rate is 
high (~ 10~^ M0yr~^) and consistent with disk-accretion 
scenarios explaining the formation of massive stars (e.g., 



cently observed the intermediate-mass young stellar object Wolfire & Cassinelli 1987 




Jijina & Adams 1996, 


Norberg 


poweriig HH288 (core mass between 6 and 30 Mq). This & Maeder 2000 




Tan & McKee 2002 


, and 


Yorke 2002|). 



outflow is quadrupolar and most likely due to two power- 
ing sources not separable by the resolution of the observa- 
tion. But in spite of being quadrupolar, the whole system 
can be explained by two outflows with shock-entrained 
material similar to HH211. 

A high-mass object: our observations of 
IRAS 05358+3543mm show an example of a mas- 
sive star-forming cluster. High-resolution Plateau de 
Bure observations resolve the single-dish features into a 
number of different outflows. So far, just a few massive 
outflows have been mapped with high spatial resolution 
(see Introduction), and we stress that the large scale 
eastern outflow is the flrst massive outflow observed 
with such a high degree of coUimation (coUimation factor 
~ 10). We estimate the accretion rate to be 10~^ M0yr~^, 
and [Beuther et al. (2002b| ) recently showed that mass 
entrainment rates of massive star-forming regions are 
usually of that order. This is high enough to overcome 
the radiation pressure of the central object and build up 



more massive stars ( 


Wolflre & Cassinelli, 1987; Jijina & 


Adamt 


, 199d; 


Norberg & Maeder, 200C ; |Tan & McKee, 


2002; 


Yorke, 2002 


). Additionally, jet models require 



disks, which are believed to be obs erved in a few massive 
objects (e.g., IRAS 20126+4104, [Cesaroni et al . 1997|, 



Zhang et al. 199^ ; IRAS 23385+6053 [MoHnari et al . 
199^). [Thus, to explain the outflow features observed 
in IRAS 05358+3543mm, shock-entrainment models are 



suflicient ( e.g., [Gueth &: Guilloteau 1999| , p^icher et al. 
2000[3nd [Cabrit et al. 1997[) , and no other formation 



mechanism is needed. It cannot be ruled out that in other 
sources different physical mechanisms are taking place, 
but these observations indicate that other massive star 
formation sites of rather chaotic appearance might be 
disentangled into simpler structures if they are observed 
at higher angular resolution. 

5. Summary 

IRAS 05358+3543mm is the flrst example of a massive 
(> 10 Mq) bipolar outflow with a high degree of col- 
Hmation on scales of 1 pc (coUimation factor ^ 10). 
High-angular-resolution observations with the Plateau de 
Bure Interferometer resolve the single-dish observations 
into at least three different outflows. Two of them form 
a quadrupolar system, most likely emanating from adja- 
cent protostars within the same mm core. Our data resolve 



CH3OH looks very similar, but the higher resolution PdBI 
observations show varying morphologies traced by the dif- 
ferent lines. While SiO is observed mainly in the outflows, 
H^^CO"'' traces core-like structures, which do not coincide 
exactly with the dust cores because of varying relative 
H^^CO+ abundances. Finally, CH3OH can be decomposed 
into a core tracing component at the line center and wing 
emission tracing the outflows. 

Ratio maps between CO 6-5, 2-1 and 1-0 reveal lo- 
cal temperature gradients. At the tip of the high-velocity 
outflow (B), we flnd a temperature increase (> 80 K) 
caused by shock-interaction of the outflow with the ambi- 
ent medium. Additionally, the southern lobe of the third 
outflow (C) is much warmer (again > 80 K) than the rest 
of the outflow, possibly due to UV heating of a close by 
and more evolved cluster. Contrasting to these tempera- 
ture increases, the possible center H^'^C0''"[3] of outflow 
(C) is cold, below 15 K. 
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